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Recent experiments have revealed spectacular transport properties of conceptually simple 
semimetals. For example, normal semimetals (e.g., WTe2) ® have started a new trend to 
realize a large magnetoresistance, which is the change of electrical resistance by an external 
magnetic field. Weyl semimetal (WSM)l^is a topological semimetal with massless relativistic 
electrons as the three-dimensional analogue of graphene ^ and promises exotic transport 
properties and surface states which are different from those of the famous topological 
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insulators (TIs)l31l In this letter, we choose to utilize NhP in magneto-transport experiments 
because its hand structure is on assembly of a WSM 013 and a normal semimetal. Such 
a combination in NbP indeed leads to the observation of remarkable transport properties, 
an extremely large magnetoresistance of 850,000% at 1.85 K (250% at room temperature) 
in a magnetic field of 9 T without any signs of saturation, and ultrahigh carrier mobility 
of 5x10® cm^ s~^ accompanied by strong Shubnikov-de Hass (SdH) oscillations. NbP 
presents a unique example to consequent design the functionality of materials by combining 
the topological and conventional phases. 

A Weyl semimetal (WSM) is a three-dimensional (3D) analog of graphene, in whieh the 
eonduetion and valenee bands eross near the Fermi energy. The band crossing point, the so-called 
Weyl point, acts as a magnetic monopole (a singular point of Berry curvature) in momentum space 
and always comes in a pairs. Unusual transport properties and surface states such as Fermi arcs 
are predicted, stimulating strong interest in realizing the WSM state in real materiald^Einil. If the 
time-reversal and inversion symmetries are respected, a pair of Weyl points is degenerate in en¬ 
ergy, forming another topological phase called Dirac semimetal E^Ell. WSMs and Dirac semimetals 
usually exhibit very high mobilities, as observed in transport experiments (such as Cd 3 As 2 )f^^^. 
Generally, semimetals are new platforms to realise a huge magnetoresistance (MR) I 13111, an effect 
that has been pursued intensively in emerging materials in recent years, because of its significant 
application in state-of-the-art information technologies Electrical transport in a semimetal usu¬ 
ally consists of two types of carriers (electrons and holes), leading to large MR when a magnetic 
field is applied with an electron-hole resonance Hill. In a simple Hall effect setup, the transverse 
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current carried by a particular type of carrier may be nonzero, although no net transverse cur¬ 
rent flows when the currents carried by the electrons and holes compensate for each other. These 
nonzero transverse currents will experience a Lorentz force caused by the magnetic field in the 
inverse-longitudinal direction. Such a back flow of carriers eventually increases the apparent lon¬ 
gitudinal resistance, resulting in dramatic MR that is much stronger than that in normal metals 
and semiconductors. Thus, it is crucial to obtain high purity samples to realize a balance between 
electrons and holes and a high carrier mobility (/i) as well, both of which will enhance the MR 
effect. 

Elemental Bi and WTe 2 (Ref. exhibit large MR as typical examples of semimetals, in 
which electron and hole pockets coexist on the Fermi surface (Fig. la). There is a special type 
of semimetal whose conduction-band bottom and valence-band top touch the Fermi surface at the 
same point in momentum (k) space (Figs, lb and Ic). Many such semimetals exhibit a high carrier 
mobility and relatively large MR with a linear dependence on the magnetic field, such as zero-gap 
topological-insulator silver chalcogenides I231I1 and Heusler compounds the Dirac semimetal 
CdgAsa (/i = 9 X 10® cm^ s'^ at 5 K, MR= 1500% at 1.5 K and 14.5 T) ElEZlill, and the 
WSM TaAs (/i = 5 x 10® cm^ s"^ at 2 K, MR= 5.4 x 10® % at 10 K and 9 T)l^. The high 
mobility may originate from the linear or nearly linear energy dispersion. The unsaturated linear 
MR is interpreted as a classical effect due to the strong inhomogeneity in the carrier density or 
as a quantum effect by the linear energy dispersion at the band touching point 1^. The semimetal 
NbP combines the main features of the WTe 2 -type (showing extremely large MR) and Cd 3 As 2 -type 
(showing ultrahigh mobility) semimetals in the band structure (Fig. Id), exhibiting hole pockets 
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from normal quadratic bands and electron pockets from linear Weyl bands. As we will see, NbP 
exhibits an ultrahigh earrier mobility eomparable to that of Cd 3 As 2 and an extremely large MR 
surpassing that of WTe 2 . 


The single erystal of NbP used for the present study and the respeetive X-ray diffraetion 
patterns are shown in Fig. 2a. The erystal strueture of NbP is non-eentrosymmetrie spaee group 
lAimd (Fig. 2b). No indieation of twinning was found in the diffraetion experiments. Both 
atom types have the same eoordination number of six, and the same eoordination environment in 
form of a trigonal prism. A detailed overview of the struetural eharaeterization is presented in 
the supplementary information. A measurement of the temperature dependenee of the resistivity, 
Pxx{T), is a simple way to identify the eleetronic states of a material. On the basis of our high- 
quality single erystals of NbP grown via ehemieal vapor transport reaetions, pxx{T) is measured 
under various transverse magnetie fields ranging from 0 to 9 T, as shown in Fig. 2e. At zero field, 
we observe metallie behavior with p^x (300 K) = 73 pVL cm and a residual resistivity pxx{‘^ K) = 
0.63 pVt em. This results in a residual resistivity ratio [pa;3;(300 K)/pxx{‘^ K)] = 115, whieh is 
direetly related to the metallieity and quality of the erystal. Compared to other similar materials 
at low temperature (2 K), NbP exhibits a resistivity that is about 30 times lower than that of WTe 2 
(Ref. [T]) but 30 times higher than that of Cd 3 As 2 (Ref. 15 1 . NbP does not beeome supereondueting 
for temperatures above 0.10 K. After applying magnetie fields, we observe a remarkable ehange in 
the resistivity. Pxx{T) ehanges from positive slope (metallic) to negative slope (semieondueting) 
from a very small field of 0.1 T and beeomes eompletely semieondueting at a field of 2 T. This may 


be due to the opening of a gap at the Weyl point. In general, a eonventional semimetal does not 
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display such behavior, whereas some small-gap or gapless semimetals (e.g., WTeJ^^and CdsAsJ^^ 
exhibit a similar trend, usually at very high fields and low temperature. Another important fact 
observed in the present material is that pxx{T) also increases dramatically due to the application 
of magnetic fields at room temperature (300 K). 

We now focus on the MR measurement in NbR The MR is commonly calculated as the ratio 
of the change in the resistivity due to the applied magnetic field {H), [{p{H) — p(0))/p(0)] x 
100%. Figure 2d displays the MR measured in transverse magnetic fields up to 9 T at different 
temperatures. At low temperatures, we find that NbP exhibits an extremely large MR = 8.5x10®% 
at 1.85 K in a field of 9 T. This MR is five times larger than that measured for the same field in 
WTeJi^ and nearly two times as large as than that of TaAs 1^, another WSM predicted in the same 
family as NbP 121121. Upon increasing temperature, the MR of NbP remains almost unchanged up to 
20 K and then starts to decrease at higher temperatures. We note that the MR is still as high as 250 
% in 9 T at room temperature (inset of Fig. 2d). We have also measured the angular dependence of 
the MR on the direction of magnetic field. Figure 2f shows the MR observed at different angles {0) 
between the field and current directions. The MR decreases slightly from 8.5 x 10® % at 0 = 90° 
(transverse) to 2.5 x 10®% at 6^ = 0° (longitudinal). Thus, the MR varies only by a factor of 3.4, 
implying the relatively isotropic nature of the material, i.e., compared to the layered semimetal 
WTea. 


A large MR is usually associated with a high mobility. The carrier mobility and concentration 
are two important parameters of a material, which can be derived from the Hall coefficient. We 
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have performed Hall effeet measurements in both temperature-sweep and field-sweep modes to 
improve the aeeuraey of our data. The field dependenee of the Hall resistivity p^y(H) exhibits 
a linear eharaeteristie at high fields (see Supplementary Information). However, the nonlinear 
behavior in low fields indieates the involvement of more than one type of eharge earriers in the 
transport properties. As seen from the inset of Fig. 2e, NbP exhibits a negative Hall eoeffieient, 
Rh{T), up to 125 K whieh ehanges sign for temperatures above 125 K. For the sake of simplieity, 
we use the single-earrier Drude band model, ne,h{T) = l/[eRH{T)], to ealeulate the earner density 
and Pe,h{T) = Rh{T) / pa^x{T) in order to estimate the mobility, where ne{nh) and Pe{^J^h) are 
eharge density and mobility of the eleetron (hole), respeetively. We use the slope of pxy{H) at high 
fields to ealeulate the Hall eoeffieient (Fig. 2e). The eleetron earrier eoneentration, Ue, is found to 
be 1.5 X 10^® em ^ at 2 K and inereases slowly with temperature, exhibiting a semimetal-like or 
very small gap-like behavior. The mobility plays a major role in the eharge transport in a material 
and eonsequently determines the effieieney of various deviees. Here, NbP exhibits an ultrahigh 
mobility of 5 x 10®om^V“^s“^ at 2 K. This value is elose to that of Cd 3 As 2 whieh is reported 
to exhibit a mobility of 9 x 10® em^ V“^s“^, one order magnitude higher than that of TaAs at 2 K. 
Furthermore, it has been shown that the mobility in Cd 3 As 2 seales with the residual resistivity. 
Henee, it ean be proposed that materials with a low residual resistivity exhibit high mobilities, 
with the present ease being a good example. 

The low field measurements (Fig. 2d and 2f) do not show any saturation of the MR up 
to 9 T in the entire temperature range. In order to pursue this finding to even higher fields, we 
have performed transverse MR measurements up to 30 T in do magnetie fields as shown in Fig. 
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3a. The MR increases up to 3.6 x 10®% for a field of 30 T at 1.3 K and still shows no tendency 
of saturation. We have further corroborated this trend of a non-saturating magneto-resistance by 
performing experiments in pulsed magnetic fields up to 62 T at 1.5 K (inset of Fig 3 a). The MR 
reaches a value of 8.1 x 10® without any indication for saturation. The MR continues increasing 
up to 62T the maximum field reached in the pulsed magnetic field experiments. At 62T we find 
an MR of 8.1 x 10®. As already indicated in the 9 T MR data depicted in Fig. 2d and 2f, SdH 
quantum oscillations appear for T < 30 K. At the lowest temperature (1.85 K) the oscillations 
start for fields as low as 1 T. Since SdH oscillations only appear when the energy spacing between 
two Landau levels is larger than their broadening due to disorder, we can estimate a lower limit to 
the quantum mobility of the carriers involved jjiq > 10® cm^V~^s“^, in agreement with the large 
electron (transport) mobility extracted from the Hall effect. 

The physical parameters of the charge carriers are derived from the Fermi surface that is 
measured by the oscillations observed in the transport properties. Both p^x and pxy measured up 
to 9 T exhibit very clear SdH oscillations starting from 1 T. This indicates a very low effective 
mass, resulting in a high mobility. To obtain the amplitude of the SdH oscillations, Apxx from pxx, 
we subtracted a smooth background. The results are plotted as a function of l/(poH) at various 
temperatures in Fig. 3b. As expected, the oscillations are periodic in 1/(pqH). They occur due 
to the quantization of energy levels which directly gives the effective mass of the charge carriers. 
A single oscillation frequency F = 32 T is identified from the data shown in Fig. 3b, which 
corresponds to 1/{pqAH) = 0.03127 T“^. We calculated a corresponding cross-sectional area 
of the Fermi surface Ap = 0.003 A from the Onsager relationship F = where 
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$0 is the single magnetic flux quantum. Supposing a circular cross section, a very small Fermi 
momentum kp = 0.0312 A ^ is obtained by assuming a circular cross section. The cyclotron 
effective mass of the carriers is determined by fitting the temperature dependence of amplitude of 
the oscillations to the Lifshitz-Kosevich formula, 


^PxxjT, B) _ -2-K'^kBTo/i3 /(5 

Pxx{^) sinh{2'K‘^kBT /13) 

where f3 = ehB/2Tim* and Tp, = h/An'^rkB are the fitting parameters, which directly 
result in the effective mass m* and quantum lifetime r of the charge carriers. The value of m* 
and r calculated from the temperature dependent SdH oscillations in 8.2 T are 0.076 mo and 
1.5 X 10“^^ s, respectively, where tuq is the bare mass of the electron. These values are comparable 
to that reported for Cd 3 As 2 We also find a very large Fermi velocity vp = hkp/{2Trm*) = 
0.48 X10® ms h The large Fermi velocity and low effective mass are responsible for the observed 
ultra-high mobility in NbP. 

The SdH oscillations at different temperatures obtained from the 30 T dc magnetic field 
measurements are plotted in Fig. 3c. At the lowest temperature, marked by arrows, the SdH 
maxima at 8 T (0.125 T-1) and 11 T(0.091 T-1) start to split into two distinct peaks, and the 
maximum at 16 T (0.0625) develops into four peaks. This can be assigned to the lifting of the 
spin-degeneracy and the degeneracy of the dual Weyl point. 

To further understand the transport properties of NbP, we have performed ab-initio band- 
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structure calculations. NbP crystallizes in a body-centered-tetragonal lattice with the nonsymmor- 
phic space group lAimd. The lack of inversion symmetry of the lattice leads to the lifting of spin 
degeneracy in the band structure (Fig. 4a). Near the Fermi energy, twelve pairs of Weyl points lie 
aside the central planes in the Brillouin zone (Fig. 4b), consistent with recent calculationd^ini. An 
important feature in the band structure is that the Fermi energy crosses the quadratic-type valence 
bands and also the linear Weyl-type conduction bands, leading to hole and electron pockets on the 
Fermi surface (Fig. 4c), respectively. Four equivalent hole pockets appear around the Z point. 
Four larger and eight smaller electron pockets exist near the S and N points, respectively, because 
the four pairs of Weyl points near S are lower in energy than the other eight pairs near N. The 
distribution of the carrier pockets is relatively isotropic compared to the layered material WTe 2 , 
for example, explaining the weak anisotropy of the observed MR. 

As we have demonstrated here, NbP is an exotic semimetal with interesting transport prop¬ 
erties. As seen in the band structure, it combines the electronic structures of a normal semimetal 
and a Weyl semimetal together. Similar to normal semimetals such as WTe 2 and Bi, NbP exhibits 
both electron and hole pockets at different positions in the Brillouin zone; however unlike WTe 2 
and Bi, its electron pockets are relevant to the linear Weyl bands, from which the high mobility 
of the high-quality samples may originate. In contrast to Weyl semimetals (e.g., TaAs) and Dirac 
semimetals (e.g., Cd 3 As 2 ), in which the Fermi energy may cross only one type of bands (electron 
or hole), NbP naturally hosts both types of carriers and consequently exhibits a huge MR in an 
electron-hole resonance situation. 
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Methods High-quality single crystals of NbP were grown via a chemical vapor transport reaction 
using iodine as a transport agent. Initially, polycrystalline powder of NbP was synthesized by a 
direct reaction of niobium (Chempur 99.9%) and red phosphorus (Heraeus 99.999%) kept in an 
evacuated fused silica tube for 48 hours at 800 °C. Starting from this microcrystalline powder, 
the single-crystals of NbP were synthesized by chemical vapor transport in a temperature gradient 
starting from 850 °C (source) to 950 °C (sink) and a transport agent with a concentration of 13,5 
mg/cm^ iodine (Alfa Aesar 99,998%) 1^. The orientation and crystal structure of the present single 
crystal were investigated using the diffraction data sets collected on a Rigaku AFC7 diffractometer 
equipped with a Saturn 724-1- CCD detector (monochromatic MoTCq, radiation, A = 0.71073 A). 
Structure refinement was performed by full-matrix least-squares on F within the program package 
WinCSD. 

The transport measurements were performed in different physical property measurement sys¬ 
tems (PPMS, Quantum Design, ACT option, home build adiabatic demagnetization stage). The 
30 T static magnetic-field measurements were performed at the High Field Magnet Laboratory 
HFML-RU/FOM in Nijmegen, and the pulsed magnetic field experiments were carried out at the 
Dresden High Magnetic Field Laboratory HLD-HZDR; both laboratories are member of the Eu¬ 
ropean Magnetic Field Laboratory (EMFL). The ab-initio calculations were performed within the 
framework of density-functional theory (DFT) with the generalized gradient approximation 1^. We 
employed the Vienna ab — initio simulation package with a plane-wave basis The core elec¬ 
trons were represented by the projector-augmented-wave potential. For calculation of the Fermi 
surface, we have used the ab — initio tight-binding method based on maximally localized Wannier 
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functions (MLWFs)l^. 
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Figure 1 Band structure for different semimetals. Schematic illustration of different types 
of semimetals and representative materials, a. Normal semimetal with the coexistence of 
electron and hole pockets, b. Semimetal with quadratic conduction and valence bands 
touching at the same momentum point, c, Weyl semimetal and Dirac semimetal, d. 
Semimetal with one hole pocket from the normal quadratic band and one electron pocket 
from the linear Weyl semimetal band for NbP. The valence and conduction bands are 
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indicated by bluish and greenish shading, respectively. The Fermi energy is marked by 


horizontal lines. 



Figure 2 Crystal structure, magnetoresistance and mobility, a, Orientation of the mea¬ 
sured single crystal NbP with the respective X-ray diffraction patterns, b, Crystal structure 
of NbP in a body-centered-tetragonal lattice, c, Temperature dependence of the resistiv¬ 
ity measured at different transverse magnetic fields displayed next to the corresponding 
curve. The inset of c shows the temperature dependence of resistance measured up to 
0.1 K in zero field, d, Transverse magnetoresistance measured at different temperatures 
with field up to 9 T. The inset shows the magnetoresistance at higher temperatures, e, 
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Temperature dependencies of the mobility (left ordinate) and carrier density (right ordi¬ 
nate). The inset shows the evolution of the Hall coefficient with temperature. The tem¬ 
perature regimes where electrons and holes act as main charge carriers are marked with 
bluish and reddish colors, respectively, f, Magnetoresistance measured at different an¬ 
gles, Q, between the current and the magnetic field as shown schematically in the inset. 





Figure 3 High-field magnetoresistance and SdH oscillation, a. Transverse MR measured 
up to 30 T static magnetic field at different temperatures. The inset of a shows MR mea¬ 
sured up to 62 T in a pulsed magnetic field, b, SdH oscillations after subtracting the 
background from the 9 T measurements. The inset of b shows the temperature de¬ 
pendence of the relative amplitude of for the SdH oscillation at 8.2 T. The solid line is 
a fit to the Lifshitz-Kosevich formula, c, SdH oscillations after subtracting the background 
from the 30 T measurements. 
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Figure 4 Structure and electronic properties of NbP. a, Ab — initio band structure. The 
Fermi energy is shifted to zero. The valence and conduction bands are indicated by the 
red and blue lines, respectively, b, The first Brillouin zone in momentum space for the 
primitive unit cell. Green and red dots represent the Weyl points with opposite chiral¬ 
ity. c, Fermi surfaces. The blue and red surfaces represent electron and hole pockets, 
respectively. 
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